Dynamic Stark effect and forbidden-transition spectral lineshapes 
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We report on an experimental and theoretical study of the dynamic (ac) Stark effect on a forbidden 
transition. A general framework for parameterizing and describing off-resonant ac-Stark shifts is 
presented. A model is developed to calculate spectral line shapes resulting from resonant excitation 
of atoms in an intense standing light-wave in the presence of off-resonant ac-Stark shifts. The 
model is used in the analysis and interpretation of a measurement of the ac-Stark shifts of the 
static-electric-field- induced 6s 2 1 So — > 5d6s 3 Di transition at 408 nm in atomic Yb. The results 
are in agreement with estimates of the ac-Stark shift of the transition under the assumption that 
the shift is dominated by that of the 6s 2 1 So ground state. A detailed description of the experiment 
and analysis is presented. A bi-product of this work is an independent determination (from the 
saturation behavior of the 408-nm transition) of the Stark transition polarizability, which is found 
to be in agreement with our earlier measurement. This work is part of the ongoing effort aimed at 
a precision measurement of atomic parity- violation effects in Yb. 

PACS numbers: 32.60.+i,32.70.-n,32.80.Ys 
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I. INTRODUCTION 

Resonant excitation of an atomic transition in an in- 
tense standing-wave light field can lead to surprisingly 
complex and non-trivial spectroscopy. Here we give a de- 
tailed experimental and theoretical analysis of the spec- 
tral line shapes resulting from off-resonant ac-Stark shifts 
when atoms in an atomic beam are resonantly excited in 
such a field. 

Spatial variation of the intense standing-wave light 
field leads to position- and velocity-dependent ac-Stark 
shifts. This results in asymmetric line shapes exhibiting 
sub-Doppler features. The attention to this problem was 
first drawn by Wieman et al. p] who studied the for- 
bidden 6S— >7S magnetic-dipole (Ml) and Stark-induced 
transition in cesium (Cs) in the context of their parity- 
violation experiment (see Ref. and references therein). 

In this work, we describe our investigation of this ef- 
fect on the 6s 2 1 So — > 5d6s 3 Di transition in atomic yt- 
terbium (Yb). While conceptually similar to the exper- 
iment done in Cs, the atomic systems and experimental 
arrangements are sufficiently different so we are able to 
study the effects in a different regime. One of the main 
differences in the two works is that the magnitude of the 
ac-Stark shifts are significantly larger in Yb than they 
are in the Cs transition studied by Wieman et al., while 
the natural line width of the forbidden 1 So — » 3 Di tran- 
sition is an order of magnitude smaller. As a result, we 
work in the regime where the ac-Stark shifts often sig- 
nificantly exceed the natural line width of the transition 
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|28j . Another key feature of this work is that the Yb 
transition studied here is between a J = ground state 
and a J' = 1 excited state while the Cs transition is 
J = 1/2 — ► J' = 1/2. Thus, in general, different tensor 
components of the ac-Stark polarizability play a role in 
these two cases. Finally, we are able to study the effect 
of the ac-Stark shifts in a regime where saturation effects 
become important. 

Our interest in this problem is related to the ongoing 
work in our laboratory on measuring parity nonconser- 
vation (PNC) in atomic Yb [j. |j. LJ [(J . Detailed under- 
standing of the line shapes with intense standing-wave 
excitation, as well as the magnitude of the ac-Stark shifts 
of the transition, are of crucial importance for an accu- 
rate measurement of parity violation in this transition. 

In this paper, we develop a theoretical treatment of 
the line shape that generalizes and extends the approach 
of Wieman et al. [ij, and is appropriate for the inter- 
pretation of our experiments with ytterbium. We also 
present a measurement of the ac-Stark shifts for the 
6s 2 1 So — > 5d6s 3 Di transition and discuss our theoret- 
ical modeling of the line shape. 



II. THEORY 

A. Tensor structure of the ac-Stark polarizabilities 

The general tensor structure of ac-Stark polarizabil- 
ities has been discussed by various authors (see e.g. 
Hll H IIS El) In these works the ac-Stark shifts are 
discussed in terms of effective static electric and magnetic 
fields. In this section, we present a parameterization of 
the ac-Stark shifts in terms of irreducible spherical com- 
ponents in a manner analogous to what is commonly done 
for dc-Stark shifts. 
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Ignoring possible time-reversal-invariance- violation ef- 
fects, there is no permanent atomic electric-dipole mo- 
ment and, consequently, there is no first-order Stark shift. 
Using time-dependent second-order perturbation theory 
and averaging the energy shift over the period of the light 
oscillation gives an energy shift for state \j, J, M) due 
to non-resonant light of |12| 

A£( 7 ,J,M)= (£o) ' 4 (£5)j J2 {l,J,M\d % W,J',M')x 

7',J',M' 

( 1 ',J',M'\d*\ 7 ,J,M)x 
1 

£(7, J,M) ~£(Y, J',M>) + huj 
1 

+ £(7, J, M) - s(y, J', Wj - 

(i) 

where the oscillating field (with angular frequency lu) is 
given by e{t) = e cos(oii), d is the dipole operator, and 
£(7, J,M) is the unperturbed energy of state I7, J, M). 
The sum is taken over all of the opposite-parity atomic 
energy eigenstates, including the continuum. 

Without referring to specific intermediate energy lev- 
els, we may say that the energy shift of a given atomic 
state, being a scalar quantity, must result from a contrac- 
tion of irreducible tensors of the same rank describing the 
light and the atom. The light tensors are bi-linear in the 
components of the light electric-field amplitude vector e 
and its complex conjugate e*. The three irreducible ten- 
sor components built from e and e* arc 
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second-rank tensor. (4) 



These three tensor components correspond to the inten- 
sity, orientation, and alignment of the light field, respec- 
tively. The vector part of the light tensor, Eq. ®, can 
be also written in terms of the dual "circular-intensity" 
vector: 



1/ * 



1 



£i£~ 



where 



(5) 



(6) 



and eijk is the totally anti-symmetric tensor. Note that 
for the circular intensity vector to be nonzero, the light 
polarization vector components cannot all be of the same 
phase. Consequently, the circular-intensity vector van- 
ishes in the limit of a linearly polarized field. As an 
example of a non-vanishing circular-intensity vector, for 
left-circular ly polarized light propagating along z, 



e >+)=-: 
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(7) 



-i\e] 2 z. 



(8) 



To obtain a shift of an atomic state, three atomic ten- 
sor operators of ranks K = 0, 1, and 2, respectively, 
must be contracted with the light tensor components 
given in Eqs. 101-101), so the ac-Stark shift operator can 
be written as 



H ac = 



+ £i£i 



£;£, 



)T 



(i) 



T 



(2) 



(9) 



Generally, the eigenstates of the Hamiltonian are su- 
perpositions of the magnetic sublevels. However, if mix- 
ing of the different magnetic sublevels is negligible, for 
example if the magnetic sublevels are split by a suffi- 
ciently strong dc magnetic field, we can find the shifts of 
the magnetic sublevels by simply evaluating the diagonal 
elements of the above Hamiltonian. In this case, only 
the tensor operators with the tensor index q = con- 
tribute to the matrix elements. These contract with the 
q = components of the electric-field tensors. By the 
Wigner-Eckart theorem, the atomic matrix elements are 
proportional to the matrix elements of operators of the 
corresponding rank and index q built from the compo- 
nents of the total angular momentum J. This results in 
an ac-Stark shift for the I7, J, M) state of 



= Co|£| 
+ C 2 ( M 2 
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(10) 



where the constants C n are determined by the specifics 
of the atomic energy levels. Here we used the explicit 
form of the tensors of rank up to 2 with q = (see, for 
example, Ref. 01, Section 3.2.2) expressed in Cartesian 
components, for example, that the q — component of 
T^isTj?. 

We choose to define the constants C„ so that the scalar 
and the second-rank tensor parts are consistent with the 
usual definitions for static electric polarizabilities: 



A cstatic 



2 



02 ( ZM 2 - J(J+1 ) ^ 3E 2 
2 V J(2J-1) 



E* 



(11) 



Here E is the static electric field. The second-rank tensor 
polarizability is defined so that the tensor shift averaged 
over M is zero and so that, if the electric field is applied 
along z, then — a^E 2 /2 is the tensor shift of the stretched 
states M — ± J. We adopt a similar convention for the 
vector term in the ac-Stark case: — ai\e] 2 /2 is the vector 
shift of the M — J state for left circularly polarized light 
propagating along z. Finally, we have a general expres- 
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This analysis allows one to define the scalar, vector, 
and tensor ac-Stark polarizabilities, but, we again point 
out that Eq. (|12fl is only applicable if the mixing of the 
magnetic sublevels due to the ac-Stark perturbation is 
negligible. If linearly polarized light is applied in the 
absence of other fields, the magnetic sublevels will split 
along an axis of quantization defined by the polarization 
of the light with energy shifts given by Eq. Ijl2|l . 

Equations © and l|12fl dictate the selection rules for 
the scalar, vector, and second-rank tensor nolaxizabili- 
ties: = 

and , 
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FIG. 1: Relevant low- lying energy levels and transitions in 
atomic Yb. The forbidden transition at 408-nm is driven by 
a standing-wave light field in a power-buildup cavity. 



B. ac-Stark Shifts on a Stark-Induced Transition 

The transition studied in this work is the highly forbid- 
den 6s 2 1 So — > 5d6s 3 Di transition in atomic Yb (Fig. 
HJ. In the absence of external fields this transition oc- 
curs through a highly suppressed magnetic-dipole (Ml) 
amplitude with a value of |5j 

ifD^MjI/^So)! = 1.33(8) x 1(T 4 M o, (13) 

where fio is the Bohr magneton. The above error is re- 
duced from that quoted in Ref. p| due to a reevaluation 
of the branching ratios for the decay of the 5d6s 3 Di state 
based on the work of Ref. . 

The ongoing PNC experiment relies upon interference 
of a parity-violating transition amplitude with a parity- 
conserving El amplitude arising from the application of a 



dc electric field. Since the transition of interest is from a 
J = state to a J = 1 state, the Stark-induced amplitude 
connecting the ground state and an M sublevel of the 
5d6s 3 Di state can be written as 



A(El st ) M = i/3(Exe)_ 



(14) 



where (E x i)_ M is the — M spherical component and (3 
is the vector transition polarizability. The magnitude of 
/3 was measured to be jj| 



2.18(10) x 1(T 8 ea /(V/cm), 



(15) 



where e is the charge of the electron and ao is the Bohr 
radius. Again, the error here is reduced from that of Ref. 
in light of Ref. 0. 

For the work discussed here, the Stark-induced El am- 
plitude is of primary interest. The nature of the Stark- 
induced transition leads to a variety of interesting con- 
clusions concerning the experimental configurations re- 
quired to access the different ac-Stark-shift parameters 
a n . We begin by considering the concrete geometry em- 
ployed in this experiment with the J = — > J = 1 
transition under study. In this configuration, the light is 
linearly polarized at some angle 9 relative the dc electric 
field as shown in Fig. [5] We define our coordinate system 
so that the axis of quantization is along the electric field 
and the polarization is nominally in the x-z plane (effects 
of misalignments are discussed in Section \V\ . 

From Eq. (f 141) we see that with this geometry, the light 
excites the state 



| 3 Di,M„ =0) = -L(| 3 D!,M Z 



-1 



\ 3 D U M Z = -1)), 
(16) 



which is a state aligned along the y axis, regardless of the 
polarization angle 0. 

This state is an eigenstate of both the dc- and ac-Stark 
perturbations, as can be seen by considering their Hamil- 
tonians, which are of the form of Eq. JJjJ. Because the 
electric field vector E does not have a y component, the 
Hamiltonian 10 is a linear combination of the operator 
products that transform as JjJj with i,j — x,z. Tak- 
ing the quantization axis along y, the matrix elements 
of these products can be written in terms of the matrix 
elements of products of pairs of the raising and lower- 
ing operators. Such products will only mix states with 
AM y = 0, 2, so the M y = state is not mixed with 
My — 1,-1. The My = state remains an energy eigen- 
state of both interactions acting simultaneously and its 
energy does not depend on the angle of the polarization 
of the light relative to the dc electric field. Based on the 
Eqs. (JTJJ, ifTTl) . and ljl"2"|l . the energy shift of this state is 



A£( 3 D 1 ,M y = 0) 



dc 



+ 2 a 



dc 



E 2 



(17) 



Here we have used the superscripts ac and dc to distin- 
guish between the static and dynamic polarizabilities (we 
omit the superscripts where the meaning of a symbol is 
clear from the context). We conclude that with the cur- 
rent configuration of fields, the experiment is only sen- 
sitive to the combination of ac-Stark parameters defined 
in Eq. I|17[l . However, it is possible to observe different 
components of the ac-Stark shifts using different config- 
urations of fields. While not directly relevant for the 
measurement described here, we briefly discuss the mod- 
ifications which would be necessary to isolate the scalar, 
vector, and tensor polarizabilities. 

Perhaps the simplest modification one can envision in 
order to separate the scalar and tensor shifts is to apply a 
strong magnetic field along the z axis. This provides the 
strong leading field needed to make Eq. (|f 2f> valid. In this 
case, the M z = — 1 and M z = +1 magnetic sublevels are 
excited independently. As can be seen from Eq. (|f 2fl , the 
ac-Stark shift of the two states will depend on the angle 
of the polarization of the light with respect to the dc- 
electric field. Measuring the shift at different polarization 
angles then allows one to separate the scalar and tensor 
components. 

In order to measure the vector ac-Stark shift parame- 
ter, one must use circularly polarized light since the shift 
depends on Im(e*e y ) [see Eqs. © and ©]. Further- 
more, the light must propagate non-orthogonally to the 
axis of quantization in order to have a nonzero compo- 
nent of the polarization-intensity vector along the quan- 
tization axis. These conditions can be achieved by using 
circularly polarized light and applying the electric field 
along the axis of light propagation. With such an ar- 
rangement, only one magnetic sublcvcl, cither M y = +1 
or My = — I, is excited and the shift of the level would 
be a combination of the scalar, vector, and tensor shifts. 



C. Estimate of Ground-State Shift 

The shift of the transition frequency is due to shifts in 
both the 6s 2 1 So ground state and 5d6s 3 Di excited state. 
Since the ground state has angular momentum J = 0, 
it can only have a scalar shift. In contrast, the excited 
state can have all scalar, vector, and tensor components. 

The ac-Stark shift is resonantly enhanced for eigen- 
states with energy separations nearby the frequency of 
the light. In addition, the effect depends quadrat- 
ically on the dipole coupling of the states involved. 
Therefore, one can expect, for light resonant with 
the 6s 2 1 So — > 5d6s 3 Di transition, the shift of the 
6s 2 1 So state to be dominated by the esBp 1 ?! state, since 
it is the nearest opposite-parity state to the 5d6s 3 Distate 
and has a strong dipole coupling to the ground state. 

The shift of the ground state can be estimated us- 
ing Eq. an( l reducing the sum to include only the 
^dGs^i) state. Using the energy separation of the 
1 Pi and 3 Di states, 579 cm -1 , and the matrix element 
|(6s6p 'Pi,M z = 0|4|6s 2 %)! = 2.4(f) ea , as deter- 
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FIG. 2: A schematic of the arrangement and orientation of 
fields for the ac-Stark experiment. 



mined from the lifetime of the 6s6p 1 Pi state [l^, we 
find 

Hz 

a (408nm, ^o) « 0.28 ^ ( 18 ) 

(V/cm) 

The shift of the 5d6s 3 Di state is significantly harder to 
estimate. It is possible that there is an odd-parity eigen- 
state with energy close to twice the energy of a 408-nm 
photon, which can lead to resonantly enhanced shifts of 
the 5d6s 3 Di state. The energy spectrum in this region 
(which is below the Yb ionization limit) is very dense 
due to the excitation of 4f orbitals [fg- While the ac- 
Stark shifts of the 5d6s 3 Di state due to the known en- 
ergy eigenstates are estimated to lead to shifts which are 
significantly smaller than that of the ground state, the 
knowledge of the energy spectrum is far from complete 
in this region. This provides one of the motivations for 
determining the ac-Stark polarizabilities experimentally. 



III. EXPERIMENTAL APPARATUS AND 
PROCEDURES 

A schematic of the experimental arrangement in pre- 
sented in Fig. |2 Atoms were resonantly excited on the 
6s 2 1 So — > 5d6s 3 Di Stark- induced transition inside an 
optical power-build-up cavity. Fluorescence at 556 nm 
from the subsequent 6s6p 3 Pi — > 6s 2 1 So transition (see 
Fig. ^ was recorded as the laser frequency was scanned 
over the resonance. The power and polarization of the ex- 
citation light were varied. A numerical calculation based 
on the density-matrix formalism was used to generate 
spectral line shapes for different values of the ac-Stark 
polarizabilities. The data were fit to these line shapes in 
order to extract a combination of the ac-Stark polariz- 
abilities for the transition. 
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Vacuum and Atomic Source 



D. Power-Build Up Cavity 



A stainless-steel oven with a multi-channel nozzle was 
used to create an effusive beam of Yb atoms inside of 
a vacuum chamber with a residual background pressure 
of w 5 x 1CP 6 Torr. The oven was heated with coiled 
tantalum-wire heaters placed inside ceramic tubes (alu- 
mina AD-998; Coors Ceramic). The oven was operated 
with the front « 100 °C hotter than the rear to avoid 
clogging of the nozzle. The typical operating tempera- 
ture was 500 °C in the rear of the oven, corresponding 
to a atomic beam density of 2.5 x 10 14 atoms/cm 3 inside 
the oven. A downstream external-vane collimator was 
used to reduce the Doppler width to « 12 MHz. The 
transparency of the collimator is estimated to be w 90% 
in the forward direction. The collimator was mounted on 
a movable platform, allowing precise alignment of the an- 
gle of the collimator relative to the atomic beam during 
the experiment. The collimator was heated using tanta- 
lum wire heaters to « 350°C to prevent clogging. 



B. Electric Field 

The dc electric held was produced by applying voltage 
to a pair of plates which were the same as those used in 
Refs. HQ. The plates were made of stainless steel and 
were 0.79 cm thick, 8.9 cm along the direction of the laser 
beam, and 3.2 cm in the direction of the atomic beam. 
The edges of the plates were rounded to a 0.3 cm radius 
and electropolishcd to avoid discharge. The plates were 
oriented so as to provide a vertical electric held and were 
separated with two Delrin spacers. The separation of the 
plates was 1.016 cm. A depression was cut into the top 
electrode and an array of holes was drilled into theplate 
in order to collect the fluorescence (see Refs. 0, Q for 
further details). 

For the work described in Ref. Q , the effect of the holes 
on the magnitude of the electric held was calculated. It 
was found that the holes reduced the electric field by less 
than 1%. The exact correction was not determined due to 
uncertainties in the surface charges of the dielectric light 
guide that was placed into the depression to measure the 
fluorescence. We therefore, conservatively estimate a 1% 
error in our knowledge of the magnitude of the electric 
held. The value of the applied field was monitored with 
a precision high- voltage divider. 



The cavity was designed as a symmetric two-mirror 
resonator with mirror radii of curvature of 25 cm. The 
cavity-mirror separation was 26.01(6) cm and the light 
was coupled in to the TEMoo transverse cavity mode. 

The cavity consisted of two high-quality mirrors 
mounted on precision optical mounts. The mirror mounts 
were attached to a Super-Invar rod supported by lead 
blocks designed to damp vibrations. 

The transmission of the mirrors was measured directly 
using a Newport 841-PE power meter with an 818-SL 
head. The incident power was measured with a calibrated 
attenuator (833-SL); the attenuator was removed for the 
measurement of the power transmitted. The power meter 
is calibrated to ±2% both with and without the attenu- 
ator. The transmission of the mirrors was measured to 
be 



T = 5.08(7)st o t(14)o«K6 x 10" 



(19) 



where the first error is from the variation of the trans- 
mission measurements and the second is the error from 
the calibration. The transmission was measured for both 
mirrors and found to be the same. 

The finesse of the cavity was measured using the 
cavity-ring-down method |l7j| . Using the experimentally 
measured ring-down time 



t = 1.17(2) ijs, 
we have a cavity finesse of 
2ttc 



T ■ 



2L 



t = 4240(70), 



(20) 



(21) 



where c is the speed of light and L = 26.01(6) cm is the 
length of the cavity. 

Using the measured finesse, we calculate a value for 
the reflectivity of the mirrors of 



1-R = 7.41(12) x 10~ 4 . 



(22) 



The total mirror loss, due to absorption and scattering, 
per bounce is estimated to be 



I = 1 - R-T = 2.4(2) x 10" 



(23) 



The amount of loss limits the fraction of power which 
can be transmitted through the cavity. The fraction of 
the power coupled into the cavity that is transmitted is 
given by 0] 



C. Residual Magnetic Field 

The magnetic held from the earth was partially com- 
pensated using three sets of magnetic coils wrapped 
around the vacuum chamber. We estimate that the 
residual magnetic held was trimmed down to a level of 
w 50 mG. 



P C \7T 



0.46(13). 



T 2 



{I + T) 



(24) 
(25) 



In order to determine the amount of power coupled into 
the cavity we monitored the light coming off of the input 
coupler of the cavity. This light consists of light that is 
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not mode matched to the cavity as well as light due to 
an imperfect cancelation of the cavity leakage field with 
the light rejected from the cavity mirror. If we couple a 
fraction of light e into the cavity, Pc — eP\ n , then the 
power coming back from the input coupler is |lq 

P R = eP in l 2 (?pJ +(l-e)P in . (26) 

The first term in this equation is due to the imperfect 
cancelation of the leakage field with the light directly 
reflected off of the mirror and represents the fundamental 
limitation on the amount of light that can be coupled into 
the cavity. The second term is the light that is rejected 
due to imperfect mode matching. 

The amount of rejected light when the cavity was on 
resonance was « 60% of that off resonance. Using this 
value and the above determinations of I and T we find 

e = 0.44. (27) 

We believe this coupling efficiency was limited by our 
inability to completely compensate the astigmatism of 
the laser beam. 

Combining the above results we find that the fraction 
of incident power that is transmitted through the cavity 
is 

^ « 0.18. (28) 

Mil 

This is in reasonable agreement with the experimentally 
measured 0.16. 

During the experiment, the power transmitted through 
the power-build-up cavity was recorded with a photodi- 
ode and the photodiode voltage was calibrated with the 
same calibrated power meter used to measure the trans- 
mission of the mirror. 



E. Laser-Frequency Locking and Stabilization 

Figure [21 shows a block diagram of the laser locking 
system. A Coherent 899 TkSaphire laser was pumped 
by sj 12 W from a Spectra Physics 2080 argon-ion 
laser operating on all lines. The Ti:Sapphirc laser pro- 
duced « 1.2 W of light at 816 nm. This light was 
frequency doubled using a commercial bow-tie resonator 
with a Lithium- Triborate crystal (Laser Analytical Sys- 
tems Wavetrain cw) to produce the 408-nm light needed 
to excite the transition. The output of the frequency 
doubler was w 80 mW. 

In order to increase the bandwidth of the frequency 
lock and narrow the line width of the laser to match that 
of the power-build-up cavity, an EOM was placed inside 
the laser cavity. The EOM used was a double-crystal 
assembly from LINOS Optics (PM 25 IR). The EOM 
was driven by a homemade electronics circuit. 

The laser was locked to the power-build-up cavity using 
the fm-sideband technique [19|. The laser beam passed 




FIG. 3: Block diagram of laser frequency stabilization. 



through an EOM (ConOptics 370) to which rf power of 
4 W at 28 MHz was coupled via a resonant LC circuit. 
The laser beam was coupled on axis to the power-build- 
up cavity. The beam passed through a glass plate which 
picked off a small fraction of the light rejected from the 
input coupler of the cavity. This light was sent to an 
amplified fast photodiode (New Focus 1801, bandwidth 
125 MHz). The signal from the photodiode was amplified 
and mixed down with the local oscillator to produce the 
error signal used to control the laser frequency. 

In order to improve absolute frequency stability, the 
resonant frequency of the cavity was locked to a her- 
metically sealed, temperature-stabilized confocal Fabry- 
Perot interferometer at 816 nm (Burleigh CFT-500) with 
a free-spectral range of 150 MHz using a similar fm- 
sideband lock as that described above. The error signal 
from this locking system was sent to a high-voltage am- 
plifier which controlled the piezo-ceramic-mounted mir- 
ror of the power-build-up cavity. Thus, the stable cavity 
provided the master frequency, with the power-build-up 
cavity serving as the secondary master for the laser. The 
frequency of the laser system was scanned over the profile 
of the atomic resonance by scanning the primary master 
Fabry-Perot cavity. Feed forward was sent to both the 
power-build-up cavity and the laser to help them track 
the frequency scans. 



F. Fluorescence Detection 

The number of atoms excited to the 5d6s 3 Di state 
was monitored by observing the fluorescence of the 
6s6p 3 Pi — > 6s 2 1 So transition at 556 nm. A lucitc light 
guide was used to send the fluorescent light to a photo- 
multiplier tube (PMT, Burle 8850) outside of the vac- 
uum chamber. A narrow-band interference filter, cen- 
tered around 560 nm, with a pass band of 10 nm was 
placed on the front of the photomultiplier tube to re- 
duce the scattered light reaching the PMT. The overall 
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detection efficiency was such that ~ 0.05% of the atoms 
undergoing the 408-nm transition produced photoclec- 
trons at the PMT photocathode. The photocathode was 
operated at —1750 V. The gain of the tube was measured 
to be w 2 x 10 6 at this voltage. The PMT current was 
fed into a Stanford Research Systems current preampli- 
fier (SRS 570). The voltage output was recorded with a 
digital oscilloscope and saved to a computer. 



G. Frequency Reference 

Since the experiment relies on a detailed understand- 
ing of the spectral line shape of the transition over a 
region of « 100 MHz, a frequency reference with closely 
spaced frequency markers is needed. To this end, we 
have constructed a Fabry-Perot interferometer operating 
at 408 nm. The mirror spacing was chosen so that the 
transverse cavity modes overlap at frequency intervals of 



2NL' 



(29) 



with N = 7 |20j. The interferometer mirrors each have 
a radius of curvature of R = 50 cm, and are separated 
by 38.9 cm. This results in a spacing between the cavity 
resonances of 55.12 MHz. 



IV. CALCULATION OF LINE SHAPE 
A. Numerical Modeling 

For a J = ^ J' = 1 transition, there are gener- 
ally four magnetic sublevels which must be included in 
a calculation. However, as discussed in Section Til Bl the 
orientations of fields used in this experiment lead to ex- 
citation of a single energy eigenstate which is a superpo- 
sition of the My = ±1 sublevels. Thus, we can neglect 
the structure of the upper state and treat the excitation 
of the atoms as a two-level problem. 

The calculation was done using the density-matrix for- 
malism. The Hamiltonian consists of the usual dipolc 
coupling between the two states with an additional term 
along the diagonal which describes the energy shifts due 
to the ac-Stark effect 



H = 







de(x, t) 



de(x,t) ujq — \ote 2 (ic) 



(30) 



Here, d is the dipole matrix element coupling the states, 
a is the ac-Stark parameter, and e(x, t) = e(x)cos(wi) 
is the oscillating electric field of the standing wave. We 
neglect the small residual running wave component of the 
field inside the cavity as it is » 1000 times smaller than 
the standing wave component. The function e(x) is the 
amplitude of the oscillating electric field as a function of 
position. Since we are considering a two-level system, 
the tensor structure of the ac-Stark interaction is not 



included. In keeping with the discussion of Section III Bl 
the a which appears here is given by 

a = (31) 
a (408nm, 3 Di) + a 2 (408nm, 3 Di) - a (408 nm, %). 

The spatial variation of the electric field translates into 
a temporal variation as an atom passes through the laser 
beam. The electric-field amplitude inside the cavity is a 
standing wave with a fundamental-Gaussian profile 



_ (£l + £f) y 

e(x)=e e r o cos(27r — ), 

A 



(32) 



where ro is the electric-field radius of the Gaussian (fol- 
lowing the convention of Ref. |2l|), and A is the wave- 
length of the light. The nominal direction of the atomic 
velocity is along x. An atom with velocity 



v = v x x + v y y + v z z , 
is subjected to an electric-field amplitude of 



(33) 



e(t) = eq e r o cos 



2vr . 

— (Vo+Vy t) 



£o e 



cos 



2tt . 

— (Vo+Vy t) 



(34) 



where we have chosen t = to be when the atom is in the 
middle of the laser beam. We have neglected v z t since 
the change in z as an atom crosses the laser beam is 
small. The largest angle an atom can have in the vertical 
direction and still intersect the laser beam is 



0i 



2L 



0.014, 



(35) 



where h = 0.64 cm is the vertical dimension of the nozzle 
and L = 23 cm is the distance between the nozzle and the 
interaction region. This means that the vertical position 
of the atom changes by only « 3% of the laser beam 
diameter as the atom crosses the beam. 

A computer code was written in the programming lan- 
guage C to numerically solve the Louiville density-matrix 
equations (see for example Ref. [2^], Ch. 15). The 
rotating-wave approximation was applied, and the equa- 
tions were solved for the upper-state population using 
the semi-implicit Euler method |23j | . The program calcu- 
lated the spectral line shape by integrating the density- 
matrix equations in time for atoms with a given trans- 
verse atomic-beam velocity v y , longitudinal atomic-beam 
velocity v x , initial zq position, and initial yo position. 
The program then integrated over each of these param- 
eters with the proper distribution functions to arrive at 
the complete spectral line shape [29). 

The upper-state population was integrated from a time 
ti = — 3\Z2r /v x to tf — +3 \/2r /v x , where ro is the 
e _1 electric-field radius. The time integral was extended 
after tf to include an exponential decay with a decay 
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time given by the lifetime of the 5d6s 3 Di state. The 
integration over the z position, or vertical position of 
the atom relative to the laser beam, was performed from 
Zi = to Zf = 3\/2fo- Calculations using these inte- 
gration ranges were found to give consistent results with 
calculations with larger integration ranges. 

The integration over the initial yo position was done 
by replacing the 2 iryo / A factor in the standing wave by 
a phase. Since the intensity of the standing wave repeats 
every 7r/2, the integration in this phase was from to tt/2. 
Atoms traveling at large velocity v y see many periods 
of the standing wave as they traverse the laser beam. 
For this reason, the initial phase has little effect on the 
resulting signal. Thus for atoms with a magnitude of \v y \ 
larger than a certain maximum value, the phase was set 
to and the signal was not integrated over the phase. 
The magnitude of the maximum velocity was empirically 
chosen to be 1.5 x Tq x X/(2tt), where To is the natural 
linewidth of the transition. 

The integration over the longitudinal velocity assumed 
a distribution given for atoms escaping from a hole |24| 

2v 3 =£ 

N(v) = —e^, (36) 
v p 

where vp — y/2kT/m w 2.7 x 10 4 cm/s is the most 
probable velocity. 

The transverse velocity distribution along y used in the 
calculation was determined experimentally. The spec- 
tral line shape of the 6s 2 1 So — > 5d6s 3 Di transition was 
recorded without the power-build-up cavity. This empir- 
ically determined line shape was then used as the input 
for the calculation. 

A significant amount of effort was put into optimizing 
the program so that it was realistic to perform numerous 
calculations with different ac-Stark shift parameters and 
electric-field amplitudes. The time it took to perform a 
calculation varied depending on the values of a and s. 
For most of the calculations done in this work the time 
was 20 — 60 min. for a single line shape consisting of 
500 frequency points, running on a 3.6 GHz Pentium 4 
processor. 

B. Calculated Line Shapes 

Experimental data were taken at a variety of different 
light powers and dc-electric fields. The maximum electric 
field of the standing wave was varied between 1.5 and 
7.5 kV/cm, and the dc electric field was varied between 
5 and 15 kV/cm. Over this range of field values the line 
shape changes quite significantly. 

Figure 0| shows the calculated line shapes for four dif- 
ferent light powers (corresponding to the amplitudes of 
the light electric field £o shown in the figure), with a 
dc electric field of 5 kV/cm. The ac-Stark parameter 
assumed for these calculations is that estimated for the 
ground-state shift, a = 0.28Hz/(V/cm) 2 . Figure[5]shows 
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FIG. 4: Calculated line shapes for different laser fields for E 
= 5 kV/cm. The ac-Stark shift used in the calculation cor- 
responds to the estimated ground-state shift. The calculated 
signals have been divided by the square of the maximum light 
electric field. The curves are vertically offset for clarity. 
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FIG. 5: Calculated line shapes for different laser fields for E 
= 15 kV/cm. The ac-Stark shift used in the calculation cor- 
responds to the estimated ground-state shift. The calculated 
signals have been divided by the square of the maximum light 
electric field. The curves are vertically offset for clarity. 



the line-shape variation with power for the same ac- 
Stark-shift parameter, but a dc electric field of 15 kV/cm. 
In both figures, the calculated signals have been divided 
by the square of the maximum electric field of the light. 

While the details of the line shapes are quite compli- 
cated, the basic shape can be understood qualitatively 
by considering several different effects. The first effect 
we discuss is the saturation of the transition. Despite 
the smallness of the transition amplitude (ss 3 x 10 _4 eao 
at 15 kV/cm), the intense standing wave present at the 
highest light powers used in the experiment is sufficient 
to saturate the transition. This is clearly observed in the 
decrease of the power-normalized signal size as the power 
is increased. The effect is more pronounced for the calcu- 
lations at a dc electric field of 15 kV/cm than for those at 
5 kV/cm due to the increased Stark- induced amplitude 
at the higher electric field. 
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In addition to the overall decrease in the signal size, 
saturation also leads to hole-burning in the velocity dis- 
tribution. The width of the dip is comparable to the 
homogeneous line width, which, in this case, is the power- 
broadened width. 

While hole-burning changes the spectral line shape, it 
does not make the line shape asymmetric around zero 
detuning. To understand the asymmetry of the line shape 
we must include the effects of the ac-Stark shift. As can 
be seen from Figs. 0]and|Sl the effect of the ac-Stark shift 
is not uniform over the profile of the line, suggesting a 
dependence on the atoms' transverse velocity. 

The source of the asymmetry can be understood in 
terms of an apparent frequency modulation resulting 
from the ac-Stark shifts. As the atoms travel through 
the laser beam the resonance frequency of the transition 
shifts depending on the intensity of the local light field. 
From an atom's perspective this is equivalent to a fre- 
quency modulation of the laser light. 

The depth of the frequency modulation (i.e., the max- 
imum frequency shift) is determined by the magnitude of 
the ac-Stark shift and is given by 



(37) 
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FIG. 6: Amplitude and intensity spectra for frequency mod- 
ulated light with different depths of modulation. The ampli- 
tudes are given by the values of the appropriate Bessel func- 
tions (see, for example, Ref. |25l|'). 



With the estimated ground-state shift and light pow- 
ers used in the experiment, the depth of the modula- 
tion ranges between 0.28 and 7.4 MHz. The frequency of 
the modulation (as the atom moves between peaks and 
valleys of the standing light wave) is determined by the 
transverse velocity 



£i 2 Vy 

2^ = ^T = 2A ^ 



(38) 



where Avjj is the Doppler shift of an atom moving with 
transverse velocity v y . The factor of two in the above 
equation arises from the fact that the ac-Stark shift is 
quadratic so that the period of the shift is half that of 
the wave length of the light. 

Since the apparent frequency of the modulation is de- 
termined by the transverse velocity [Eq. (|38|l ]. the effect 
is quite different for atoms with different v y . The trans- 
verse velocities where the character of the frequency mod- 
ulation changes, depend on the size of the ac-Stark shift, 
and, as discussed below, on the homogeneous line width 
of the transition. 

The character of frequency modulation is determined 
by the relative size of the depth of the modulation and 
the frequency of the modulation. In the case where 
£/2 << f2, frequency modulation is characterized by 
the appearance in the spectrum of the light-frequency 
detuning as perceived by the moving atom of principal 
sidebands located at ±f2 of the fundamental frequency. 
The sign of the amplitude of the sidebands is opposite 
for the high- and low- frequency sidebands. As the depth 
of the modulation is increased, the sidebands increase 
and higher-order sidebands develop. Eventually, the bulk 



of the amplitude shifts away from the fundamental fre- 
quency to the extreme frequencies. This evolution is 
shown in Fig. [S] for the electric-field amplitude and in- 
tensity of the light (proportional to the square of the 
amplitude) . 

An atom with a transverse velocity v y = Avp A sees 
two "carrier" frequencies, 



v 1 = v - Av D 
v-2 = vq + Av D , 



(39) 
(40) 



corresponding to the Doppler shifts of the two counter- 
propagating waves which make up the standing wave. 
Since the resonant frequency of the transition is modu- 
lated, the detuning of the atomic resonance relative to 
the light frequency acquires sidebands. From the atom's 
perspective, this appears as sidebands spaced a distance 
2 Ave, around v\ and v^. Thus, the high-frequency side 
band of v\ overlaps with v-i and the low-frequency side 
band of overlaps with v\. However, since the high- 
frequency sideband for v\ is positive, while the low- 
frequency sideband of is negative, there is a reduction 
of the amplitude at v\ and an enhancement at vi . This 
leads to asymmetries in the resonant peaks at ±Az/d. As 
the transverse velocity is increased, the amplitude of the 
sidebands decreases and the height of the two Doppler- 
shifted resonances evens out. 

Atoms with small transverse velocity pass through the 
laser beam without traversing multiple nodes and anti- 
nodes of the standing wave. The nature of the mod- 
ulation is therefore not sinusoidal for these atoms, and 
cannot be parameterized in terms of fl. The atoms 
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experiencing this regime have transverse velocities cor- 
responding to Doppler shifts less than the transit line 
width which is « 200 kHz. Furthermore, if SI is less than 
the homogeneous line width (determined in this case by 
the power-broadened width), then the effects of the fre- 
quency modulation are smeared out. For these reasons, 
the regime of harmonic modulation where f2 is much 
smaller than £/2 is never realized in the system under 
consideration. 

Even if the atoms move exactly perpendicular to the 
light, there is spectral asymmetry due to the spatial de- 
pendence of the ac-Stark shift - the shift is different de- 
pending on the phase of the light wave the atoms see 
(determined by the y coordinate), and on the position of 
the atoms in the x-z plane. 

The resulting spectrum is the sum of all velocity 
classes, this sum leads to asymmetries near the center 
of the spectra shown in Figs. 0]and[21 

V. ANALYSIS AND RESULTS 

Line shapes of the type shown in Figs. 0] and [S] were 
generated over the appropriate range of transition am- 
plitudes and a range of ac-Stark-shift parameters. The 
calculated line shapes were used to create an interpolat- 
ing function in frequency, the ac-Stark-shift-parameter, 
and transition amplitude. 

The calculations were related to the experimental spec- 
tral line shapes using the following parameterization of 
the calculations 

C(v) = a f(e , a, ft E, 6, s{v - v c )), (41) 

where a is the amplitude of the signal, and the function 
f(so, a, (3, E, 9, s(y — v c )) is the calculated line shape 
for a maximum electric field of the optical standing wave 
of £o, ac-Stark shift polarizability a, Stark-transition po- 
larizability ft dc electric field E, polarization angle 8, 
and frequency s(v — u c ), centered around v c . The fre- 
quency term includes a multiplicative scaling factor s. 
This s parameter was introduced in order to account for 
any variation in the overall width of the data relative to 
the calculation. Such a variation might arise due to a 
misalignment of the atomic beam relative to the axis of 
the power-build-up cavity or a deviation of the velocity 
distribution from that measured without the cavity. We 
note that the parameters in Eq. I|41[) are not indepen- 
dent. Aside from scale and offset parameters, the calcu- 
lated line shape depends only on the Rabi frequency of 
the transition (Qr = (3 E e^smO) and the ac-Stark shift 
(AS = -1/2 a eg). 

Approximately twenty experimental scans were per- 
formed with a given dc-clcctric field, polarization angle, 
and light power. These ~ 20 scans were grouped together 
to form a set of data which was analyzed as described be- 
low. The dc-electric field was varied from 5 to 15 kV/cm 
between sets. Polarization angles of 90° and ±45° rela- 
tive to the dc-electric field were used. The light power 



dc Electric Field 


ac-Stark Shift 


x 2 


N 


5 kV/cm 


-0.3187(51) 


0.77 


4 


8 kV/cm 


-0.3165(31) 


0.46 


8 


10 kV/cm 


-0.3229(20) 


1.99 


16 


Pol. = 0° 


-0.3170(46) 


- 


1 


Pol. =45° 


-0.3211(31) 


0.973 


10 


Pol. = -45° 


-0.3282(33) 


4.12 


5 


12 kV/cm 


-0.3125(29) 


0.82 


2 


15 kV/cm 


-0.3039(16) 


3.21 


15 



TABLE I: Extracted values of ac-Stark shift for different dc 
electric fields and angles of light polarization with respect to 
the dc electric field. The value of the Stark transition polariz- 
ability assumed in this analysis is /3 = 2.18 eao/(V/cm). The 
units for the ac-Stark shift are Hz/(V/cm) 2 . 



was varied from 180 nW to 2.4 mW transmitted (cor- 
responding to a maximum light electric field inside the 
cavity of 2 to 7.4 kV/cm). Data were taken at a total of 
45 different field parameters. 

The maximum light electric field inside the cavity, Eq, 
was determined using the measured power transmitted 
through the cavity and the experimentally determined 
transmission coefficient of the cavity mirrors [Eq. (|19|) ]. 
The dc-electric field was measured using the voltage 
divider mentioned in Section IIII Bl The value of the 
Stark transition polarizability was fixed in the analysis 
at /3 = 2.18 x 10~ 8 ea /(V/cm) [see Eq. ljTo|l]. Effects 
arising from the uncertainty in ft as well as an indepen- 
dent evaluation of f3 from the present data are discussed 
below. 

The power-normalized scans were individually fit to 
Eq. I|41|) . Typical scans and fits for three different com- 
binations of dc electric field and laser light power are 
shown in Fig. For a given power and dc electric field 
the mean fit values of a set of scans were determined, 
and the distribution of the values was used to estimate 
the statistical error. The overall result, including sets at 
all of the dc-electric fields, polarization angles, and light 
powers used in the experiment is 

ag c ( 3 D 1 )+<( 3 Di)-«S c ( 1 So) 

= -0.3123(11) Hz/(V/cm) 2 , (42) 

where the error is determined from the statistical errors 
in each set of scans. The reduced chi squared for the 
data is 3.18 for 45 points. The large reduced chi squared 
is primarily due to inconsistency in the data taken at 
a dc electric field of 15 kV/cm with the data taken at 
lower fields. Table [I] shows the results for the different dc 
electric fields and Fig. Q3] shows the dc-field dependence 
of the ac-Stark shift parameter. 

Excluding the data taken at 15 kV/cm gives a value of 
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FIG. 7: Experimental power-normalized line shapes and fits 
to the model described in the text for three different combina- 
tions of dc and light electric field values. Typical time for each 
scan was 200 ms. For the eo = 6.8 kV/cm (2 mW of trans- 
mitted light power), E=15 kV/cm data, the peak detected 
PMT-anode current was « 1.2 [iA. 



the ac-Stark shift of 

a§ c ( 3 D 1 ) + ^ c ( 3 D 1 )-<( 1 So) 

= -0.3188(14) Hz/(V/cm) 2 , (43) 

with the reduced chi squared of 1.57 for 30 sets. Exclud- 
ing everything but the 15 kV/cm data gives 

c^( 3 D 1 ) + e£ c ( 3 D 1 )-ag c ( 1 So) 

= -0.3039(16) Hz/(V/cm) 2 , (44) 

with a reduced chi squared of 3.21 for 15 sets. 

The discrepancy is also apparent in the s parameter 
described above. The s parameter has a clear dependence 
on the dc-electric field, going from s — 0.984(3) at 5 
kV/cm to s = 0.9535(6) at 15 kV/cm. A value of the s 
parameter less than unity corresponds to an experimental 
line shape which is broader than the line shape calculated 
assuming s = 1. 

This broadening of the line shape might suggest that 
the inconsistency in the data is a result of not properly 
accounting for saturation effects. However, there is no 
clear dependence of the extracted ac-Stark shift value on 
the light power, which one would expect if the saturation 
effects were the source of the discrepancy. 
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FIG. 8: The magnitude of the ac-Stark shift parameter 
as a function of dc electric field. The value of the Stark 
transition polarizability assumed in this analysis is /3 = 
2.18 eao/(V/cm). The error bars in this plot are based on 
the spread of the extracted ac-Stark parameters for a given 
value of the dc electric field. 



Analysis Conditions 


ac-Stark Shift 


x 2 


N 


P = 2.18 x 10" 8 ea /(V/cm) 


-0.3123(11) 


3.18 


45 


Excluding 15 kV/cm data 


-0.3188(14) 


1.57 


30 


Only 15 kV/cm data 


-0.3039(16) 


3.21 


15 


P = 2.08 x 10~ 8 eao/CV/cm) 


-0.2953(11) 


7.53 


45 


Excluding 15 kV/cm data 


-0.3112(14) 


1.88 


30 


Only 15 kV/cm data 


-0.2791(14) 


1.93 


15 


P = 2.28 x 10" 8 ea /(V/cm) 


-0.3290(11) 


2.46 


45 


Excluding 15 kV/cm data 


-0.3292(14) 


1.92 


30 


Only 15 kV/cm data 


-0.3286(19) 


3.76 


15 


Misaligned Velocity Distribution; 








P = 2.18 x 10~ 8 ea /(V/cm) 


-0.3076(11) 


10.23 


45 


P = 2.08 x 10" 8 ea /(V/cm) 


-0.2930(10) 


18.46 


45 


P = 2.28 x 10~ 8 ea /(V/cm) 


-0.3259(12) 


5.25 


45 


's' Parameter fixed at unity 


-0.3384(10) 


9.16 


27 


Atomic velocity vp = 30cm/ms 


-0.2963(11) 


7.77 


45 



TABLE II: Extracted values of ac-Stark shift for different 
analysis procedures and assumptions. The units of the ac- 
Stark shift are Hz/(V/cm) 2 . For each case, the value of the 
reduced x 2 , an d the number of sets N used in the analysis are 
presented. For all cases except the one presented in the last 
line, the most probable velocity in the beam is assumed to be 
vp — 27 cm/ms. 



We analyzed the data in a number of different ways in 
order to understand this discrepancy as well as estimate 
systematic errors from calculational assumptions. The 
results of these investigations are tabulated in Tab. |n] 
and are discussed below. 

The ac-Stark shift and the value of the Stark transition 
polarizability, /3, are highly correlated in the parametriza- 
tion of the spectral line shape. The effects of the uncer- 
tainty in the value of (3 were investigated by reanalyzing 
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FIG. 9: The magnitude of the ac-Stark-shift parameter as a 
function of dc electric field for different values of the transition 
polarizability /3. 



the data with the procedure described above assuming 
different values of j3. The results of this analysis are 
shown in Tab. |H] Figure shows the dependence of 
the ac-Stark shift parameter on the dc electric field for 
the different values of the Stark transition polarizabil- 
ity. This analysis reveals a definite dependence of the 
ac-Stark shift parameter determined from the analysis 
on the value of (3. Furthermore, the discrepancy between 
the 15 kV/cm data and the lower field data was reduced 
when the data were analyzed with a value of the Stark 
transition polarizability of /3 = 2.28 x 10 -8 eao/(V/cm). 
However, the data are still inconsistent with reduced chi- 
squared of 2.46 for 45 points. This analysis suggests 
that the value of the Stark transition polarizability may 
be higher than the central value in Eq. I] 15(1. although 
consistent within the quoted uncertainty. Below, we de- 
scribe a procedure for determining the value of the Stark 
transition polarizability from these data. This analysis 
provides a determination of the uncertainty in the ac- 
Stark shift parameter resulting from the uncertainty in 
P. Based on the results presented in Tab. ITU we arrive at 
an uncertainty in the ac-Stark shift due to uncertainty in 
(3 of 5%. 

Another possible source of experimental error arises 
from uncertainties in the atomic velocity distribution 
along the axis of the power build-up cavity. This uncer- 
tainty is the motivation for introducing the s parameter 
described above and is highly correlated with it. The 
fact that the s parameter is less than unity suggests that 
the center velocity of the atomic beam is not exactly 
perpendicular to the cavity axis. To analyze the effects 
of such a misalignment, the calculation described above 
was modified to include a Doppler distribution which was 
misaligned from the axis of the power-build-up cavity by 
0.14°. This misalignment corresponds to a Doppler shift 
of the two counter-propagating light waves constituting 
the standing wave by ±2 MHz, respectively, as seen from 
the point of view of the atoms in the atomic beam moving 
with the most probable velocity. Misalignments larger 



than this lead to noticeable effects in the velocity distri- 
bution which were not present in the experimental data. 
We believe this to be an upper limit on the possible mis- 
alignment. The data were analyzed as before; the results 
are presented in Tab. [H] From these we determine that 
the uncertainty in the ac-Stark shift parameter due to 
the uncertainties in the alignment of the atomic beam 
with the cavity axis are 2%. The values of the s parame- 
ter resulting from this analysis ranged from 1.05 for E=5 
kV/cm to 1.01 for E=15 kV/cm. 

We also investigated the effect of including the s pa- 
rameter in the fit. We analyzed a fraction of the data 
(27 sets) fixing the s parameter to unity. The fits using 
this procedure were not as good as the fits where s was 
allowed to vary. In this case, the values of the ac-Stark 
shift parameter are largely determined by the wings of 
the Doppler distribution, since that is the only parame- 
ter available to broaden the line shape. The result of this 
analysis is shown in Tab. [H] 

The longitudinal velocity distribution used in the cal- 
culation is also a possible source of uncertainty. We as- 
sume a standard longitudinal velocity distribution of ef- 
fusive atoms escaping from a hole of [Eq. I|36[)]. with the 
temperature fixed by the measured temperature at the 
rear (i.e., the coldest part) of the oven. However, atomic 
beams frequently deviate from this textbook distribution 
due to preferential scattering of the slow atoms (the so- 
called Esterman effect, see, for example, a discussion in 
Ref. 26]). In addition, the oven is not at a uniform 
temperature; the temperature varies by over 100 °C from 
front to back of the oven. This variation in tempera- 
ture can modify the observed velocity distribution of the 
atomic beam. In order to investigate the effects of uncer- 
tainty in the longitudinal velocity distribution, data were 
analyzed using calculations assuming velocity distribu- 
tions corresponding to different temperatures. The ex- 
tracted values of the ac-Stark shift parameter are shown 
in Tab. [H] Assuming a maximal deviation in the central 
temperature of 100°, we arrive at an uncertainty in the 
ac-Stark shift parameter of about 5%. 

Another key assumption in the calculation is that the 
problem can be reduced to that involving just two atomic 
levels. If the dc electric field and light field are the only 
fields present and are aligned appropriately, this is a cor- 
rect assumption (see Sect. II V|) . However, if there are 
additional fields present, such as a magnetic field, this as- 
sumption may break down. With the experimental geom- 
etry used here (Fig.|2Jl the strong dc electric applied along 
z splits the M z = ±1 3 Di states from the M z = state. 
A residual field along the x- or y-axes, would lead to mix- 
ing between the M z = ±1 and M z = states. However, 
since the Zeeman shifts (At/ = g^B w 35 kHz) corre- 
sponding to the residual magnetic fields are small relative 
to the dc-Stark shifts, (Au = 1/2 a 2 E 2 w 100 - 800 kHz 
0), this mixing is suppressed due to the strong dc elec- 
tric field al ong the z-axis and the two-level approximation 
still holds 30]. However, a residual magnetic field along 
the z-axis would serve to split the degenerate M = ±1 
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FIG. 10: Value of the extracted ac-Stark shift parameter as 
a function of magnetic field in the vertical direction. The dc 
electric field for these data ranged between 5 kV/cm and 15 
kV/cm. 
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FIG. 11: Values of the extracted ac-Stark-shift parameter as 
a function of polarization angle. The dc electric field was 10 
kV/cm for these data. 



states and the two- level treatment would be incorrect. 
To investigate this possibility, we utilized the magnetic- 
field coil normally used to compensate the ambient field 
in the z-direction to apply a magnetic field parallel to the 
dc electric field. Data were then taken at different val- 
ues of magnetic fields, which were analyzed as described 
above. The results are shown in Fig. At large values 
of the magnetic field (~ 1 G), there is a dependence in 
the extracted value of the ac-Stark shift parameter as ex- 
pected. However, no dependence was observed for data 
taken at magnetic field values that might realistically be 
present in the experiment. Thus, using the errors for 
the data taken with varying magnetic field we estimate 
the uncertainty in the ac-Stark shift parameter due to 
residual magnetic-field effects to be < 3%. 

Misalignment of the optical field relative to the dc- 
electric field can also lead to a splitting of the energy 
eigenstates that are excited. With the geometry used in 
this experiment (Fig.EJ, there are two possible kinds of 
such misalignment. One is an error in the angle between 
light field and dc electric field. Since the ac-Stark shift 
does not depend on the angle of the light field relative 
to the dc field in this geometry (see Sect. Ill Bp . a mis- 
alignment of this kind does not lead to a breakdown of 
the two-level approximation, but changes the transition 
matrix element [Eq. (|14() ]. In the other kind of misalign- 
ment, the cavity axis has a component along the dc elec- 
tric field. In this case, the assumption of a two-level sys- 
tem generally breaks down. However, for the case where 
the light field is polarized along the £-axis, the dc electric 
field and the light field remain perpendicular. Thus, the 
above argument that was used to reduce the problem to 
that involving just two atomic levels is valid in this case, 
and the misalignment does not affect the extracted value 
of the ac-Stark-shift parameter. For a general polariza- 
tion, a misalignment could influence the ac-Stark shift 
measurement. We can place a limit on the uncertainties 
arising from such possible misalignment by comparing 
the extracted ac-Stark-shift parameters for data taken 



Parameter Uncertainty Effect on a 



P 


5% 


5% 


Transit Velocity 


5% 


5% 


Calculational Approximations 


n/a 


5% 


dc Electric Field 


1% 


1% 


Polarization alignment 


w 2° 


2% 


Residual Magnetic Field 


< 50 mG 


3% 


Mirror Transmission 


3% 


3% 


Power 


4% 


4% 


Total Uncertainty 




11% 



TABLE III: Errors and the resulting uncertainty in the value 
of the ac-Stark shift parameter. 



with the angle of the light polarization relative to the dc 
electric field of ±45° with data taken at an angle of 0° . 
FigurelTTlshows the value of the ac-Stark shift parameter 
as a function of the polarization angle for data taken at 
a dc-electric field of 10 kV/cm. Using these results we 
are able to place a limit on the uncertainty in the ac- 
Stark-shift parameter due to misalignment of the light- 
polarization of 2%. 

Additional errors in the data arise from uncertainties in 
the dc-electric-field and light-field values. These errors, 
along with those discussed above are listed in Tab. HU 
The uncertainty in polarization alignment was estimated 
based on the relative signal sizes at the different polar- 
ization angles. The uncertainty of » 2° is comparable 
to the measured degree of ellipticity in the transmitted 
light due to a birefringence of the cavity mirrors. 

Including all of the systematic errors listed in Tab. IIIII 
we arrive at a systematic error of 11% in the ac-Stark 
shift parameter, while the spread of the data is « 6%. 
We arrive at a final value of the ac-Stark transition po- 
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dc Electric Field ac-Stark Shift 



N 



Parameter 



Uncertainty Effect on j3 



5 kV/cm 

8 kV/cm 

10 kV/cm 
Pol. = 0° 
Pol. = ±45° 

12 kV/cm 

15 kV/cm 

All Data 



-0.3261(38) 2.301(18) 

-0.3196(11) 2.166(5) 

-0.3289(11) 2.235(3) 

-0.3327(9) 2.283(8) 

-0.3248(13) 2.233(4) 

-0.3094(13) 2.203(5) 

-0.3284(5) 2.237(2) 



0.9931(26) 4 

0.9724(8) 8 

0.9788(21) 1 

0.9790(6) 15 

0.9654(9) 2 

0.9622(5) 15 

0.9730(3) 45 



TABLE IV: Extracted values of ac-Stark shift, Stark tran- 
sition polarizability, and 's' parameter as determined by 
chi-squared fits to the indicated data sets. The units of 
the ac-Stark shift are Hz/(V/cm) 2 and the units of (3 are 



E 


1% 


1% 


Polarization alignment 


« 2° 


1.2% 


Mirror Transmission 


3% 


1.5% 


Power 


4% 


2% 


Other 




^5% 


Total Uncertainty 




6% 



10- 8 ea /(V/cm). 



TABLE V: Systematic errors and the resulting uncertainty in 
the value of the Stark transition polarizability. The "Other" 
errors include those from uncertainties in transit velocity, cal- 
culational approximations, and residual magnetic field. These 
contributions have been estimated indirectly from their effects 
on the ac-Stark parameter (Table lTTTI and the correlation be- 
tween a and /3 determined by modeling. 



larizability of 

«H 3 D 1 )+ar( 8 D 1 )-aS c ( 1 S ) 

= -0.312(34) Hz/(V/cm) 2 . (45) 

In an attempt to determine the value of \(3\ from the 
data, we also analyzed the results using a global fit of 
all data. In addition, the values of the Stark transition 
polarizability [Eq. I|15[)] and the frequency scaling pa- 
rameter s directly affect the extracted value of the ac- 
Stark polarizabilitics. In order to properly account for 
this correlation, the variance of the data with the calcu- 
lation was simultaneously minimized with respect to all 
three of these parameters. The data were analyzed by 
varying a, (3, and s and fitting the individual data scans 
to the amplitude, a, and the center position v c . For a 
given set of a, f3, and s the \ 2 °f the individual fits were 
summed. This gave a global "x 2 map," using which a 
minimum was found with respect to a, /3, and s. 

The value of the minimum of the global x 2 was nor ~ 
malized to unity. The errors were assigned by deter- 
mining the deviations in a, (3, and s required to give 
a x 2 — 1 + 1/-^) where N is the total number of data 
points used in all of the fits, and taking the appropriate 
projections along a, (3, and s. The results for the data 
grouped by field configurations as well as the result of a 
global fit to all of the data are shown in Tab. IIVI 

There is good agreement between the results deter- 
mined from the chi-squared map and the result deter- 
mined from the fits to the individual data sets. Again, 
there is a discrepancy between the different electric-field 
sets. In addition, the extracted value of the Stark tran- 
sition polarizability (3 is inconsistent among the different 
electric field sets. Using the spread as an estimate of 
the uncertainty, we arrive at a measurement of the Stark 
transition polarizability of 

\(3\ = 2.2At° il x KT 8 ea /V/cm, (46) 

which is in good agreement with and of comparable un- 
certainty as the previous measurement [Eq. i|15|) ]. The 



error budget for the determination of (3 is presented in 
Table [V] Considering the two measurements indepen- 
dent, we get the final value 

\(3\ = 2.19(8) x 1(T 8 e a /V/cm. (47) 



VI. CONCLUSION 

The results of the measurement of the ac-Stark shift of 
the 6s 2 % -> 5d6s 3 Di transition at 408 nm [Eq. gSJl] 
are in good agreement with estimates of the ac-Stark 
shift of the 6s 2 1 So ground state [Eq. JTSJ]. The value of 
the ac-Stark shifts of this transition should not limit our 
ability to study parity nonconserving effects: the broad- 
ening due to ac-Stark shifts will not preclude resolving 
Zeeman components of the transition as necessitated by 
the current PNC-experiment scheme, and the present re- 
search has shown that lineshapes can be understood and 
modeled (see Fig. 0. In the future high-precision PNC 
experiments it will be necessary to work with a higher- 
finesse cavity. In that case, it may be possible to reduce 
the ac-Stark shifts by operating the cavity in the confocal 
regime, and expanding the laser beam to simultaneously 
excite a large number of the degenerate transverse modes. 

An additional result of the present work is a new in- 
dependent determination of the Stark transition polar- 
izability (3. The result [Eq. is consistent with the 
earlier result [Eq. ^3 of Ref. [4( and has comparable un- 
certainty, but it does not rely on the knowledge of the 
decay branching ratios of the 5d6s 3 Di state. 
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